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ABSTRACT 

We present new deep VLT/FORS optical spectra with intermediate resolution and 
large wavelength coverage of the compact radio source and ULIRG PKS 1345+12 (4C 
12.50; z = 0.122), taken with the aim of investigating the impact of the nuclear activity 
on the circumnuclear interstellar medium (ISM). PKS 1345+12 is a powerful quasar 
(L(H/3)nlr^10 42 erg s _1 ) and is also the best studied case of an emission line outflow 
in a ULIRG. Using the density sensitive transauroral emission lines [S II]4068,4076 
and [O II]7318, 7319, 7330, 7331, we pilot a new technique to accurately model the 
electron density for cases in which it is not possible to use the traditional diagnostic 
[S II]6716/6731, namely sources with highly broadened complex emission line profiles 
and/or high (N e > 10 4 cm -3 ) electron densities. We measure electron densities of N e = 
(2.941?;^) x 10 3 cm" 3 , N e = (1.47+° f 7 ) x 10 4 cm" 3 and N e = (3.16+^?) x 10 5 cm^ 3 
for the regions emitting the narrow, broad and very broad components respectively. 
We therefore calculate a total mass outflow rate of M = 8lg yr -1 , similar t o the 
range estimated for another compact radio source, PKS 1549-79 (|Holt et al.l |2006t). We 
estimate the total mass in the warm gas outflow is M to t a i = (8+|) x 10 5 Mq with filling 
factors of e = (4.4^15) x 10 -4 and e = (1.6^05) x 10~ 7 for the regions emitting the 
broad and very broad components respectively. The total kinetic power in the warm 
outflow is -Etotai = (3. 4+^3) x 10 42 erg s _1 . Taking the black hole properties published 

bv lDasvra et all (|2006t ). we find that only a small fraction {E/L hol = (1.3 + 0.2) x 10~ 3 ) 
of the available accretion power is driving the warm outflow in PKS 1345+12, which 
is significantly less than that currently required by the majority of quasar feedback 
model s (~5-10% of I/boi), but similar to the recent suggestion of iHopkins fc Elvii 
(2010) if a two-stage feedback model is implemented (~0.5% of -Lboi)- The models also 
predict that AGN-driven outflows will eventually remove the gas from the bulge of the 
host galaxy. Our observations show that the visible warm outflow in PKS 1345+12 
is not currently capable of doing so. However, it is entirely possible that much of the 
outflow is either obscured by a dense and dusty natal cocoon and/or in cooler or hotter 
phases of the ISM. This result is important not just for studies of young (GPS/CSS) 
radio sources, but for AGN in general. 

Key words: ISM: jets and outflows - ISM: kinematics and dynamics - galaxies: 
active - galaxies: ISM - galaxies: kinematics and dynamics - galaxies: individual: PKS 
1345+12 (4C12.50) 



1 INTRODUCTION 

It has become increasingly clear that the evolution and 
growth of supermassive black holes and their host galaxy 
bulges are intricately linked, through the discovery of tight 
* E-mail: jholt@strw.leidenuniv.nl correlations between the bulge and black hole proper- 
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ties (e.g. iFerrarese fe Merrittl l2000l; iGebhardt et all l200d: 



Tremaine et al.ll2002i ; iMarconi fc Hunj|2003l ; iBenson et all 
However, whilst the observed correlations become 



better constrained, the p hysical processes resp onsible are 
still not understood (e.g. Cattaneo et al.l 120091 ') with both 



AGN feedback (e.g. ISilk fc Rees 



pernova feedback (e.g. IWild et al 



19981 ) and starburst su- 
|2010| ) providing plausible 



outflow driving mechanisms. 

A number of theoretical models have successfully re- 
produced the observed correlations by including AGN 
feedback in which the rapidly accreting black hole dissi- 
pates energy into the surrounding ISM by driving power- 
ful outflows, which h alt further accretion and star forma- 
19981 ; iFabianl Il999l ; |Pi Matteo et all 
Currently, these models re- 



tion ( e.g. [S ilk & Rccs 



2005; ISpringel et al 



2005). 



quire a significant amount of the accretion energy to be 
mechanically coupled to the ISM (~5-10% of Lboi; e.g. 



Pi Matteo et al.ll2005l; iTremaine et al.ll2002l; ISpringel et all 



20051 : iKurosawa et al.M2009l: iBooth fc Schavd 120091 ). al- 
though recent work by [Hopkins fc Elvij i 201Cf ) suggests that 
if a two-stage feedback model is implemented, the initial en- 
ergy requirement may be a factor of 10 lower. 

Over the last decade, the signatures of such mechan- 
ical AGN feedback (outflows) have been observed in Com- 
pact Steep Spectrum (CSS; Dradio < 15kpc) and GigaHertz- 
Peaked Spectrum (GPS; -D ra dio < lkpc) radio sources, 
both in the war m (e.g. optical emission line outflows; e.g. 
iHolt et all 20081) and cold (e.g. HI absorption line outflows; 
e.g. iMorganti et al.l 120051 ') gas. In the infra-red, emission 
line outflows have al so recently been discovered in ULIRGs 
l[Spoon fc Holtll2009l ). 

CSS and GPS radio sources are ideal objects for study- 
ing the effects of AGN feedback as they are b elieved to be 
young, recently triggered radio-loud AGN (see IO'Dealll998l 
for an extensive review of these sources). The similarities in 
scale of the AGN, compact radio source and the dense cir- 
cumnuclcar ISM mean that interactions will be strong and 
readily detectible. In 2006, Holt et al. performed a detailed 
study of the outflow in the southern compact flat-spectrum 
radio source PKS 1549-79, deriving a mass outflow rate of 
0.12 < M < 12 M yr -1 and a kinetic power of 5.1xl0 40 

< E < 5.1 x 10 42 erg s -1 , which accounts for only a small 
fraction of the available accretion energy in this source (10 -6 

< E/L-Edd ~ 10 -4 ). However, due to the large uncertainties 
in the measured density, it was impossible to derive these 
crucial parameters more accurately. 

The young radio galaxy P KS 1345+12 (GPS; D ra dio 
~350pc; IStanghellini et al.lll993n is a prime object for AGN 
feedback studies as it contains all of the signatures of a 
recently triggered powerful AGN currently shedding its 
natal cocoon. Fast outflows (up to ~2000km s -1 ) are de- 
tected in both t he warm and c old circumnuclear ga s (e.g. 
IHolt et al.ll2003l : H03 hereafter. iMorganti et ai1l2005l ). The 
host galaxy is highly disturbed; the double nucleus (sepa- 



1 L(H/3)nlr~1 42 erg s -1 ; c.f. the SDSS DR7 quasar catalogue 
(Schneider 2010) which includes H/3 luminosities split into the 
NLR and BLR components. The NLR H/3 luminosity of PKS 
1345+12 (we no not detect the BLR) is well above the mean value 
of NLR H/3 luminosities in quasars, making it a very powerful 
AGN. This is also evident in its large [O III] luminosity. 



rated by ~1.8 arcsec or 4.3 kpc) is embedded in an extended, 
asymmetric halo with clear evidence for ex tended, distorted 
morphology (e.g. ISmith fc Heckman|[l989h . suggesting PKS 
1345+12 has been involved in a major merger in its recent 
past. In addition, prodigious star formation activity is ob- 
served. PKS 13 45+12 has a substa ntial far-IR excess {Lfir 
= 1.7xl0 12 Lm. lEvans et al.| [T999), qualifying this source as 
an ultra-luminous IR galaxy (ULIRG). Furthermore, signif- 
icant young stellar populations hav e been detected in the 
halo, both in the diffuse emission dTadhunter et al l 120051 ) 
and in a number of young Super Star Clusters (SSCs; 
iRodriguez Zaurin et al.ll2007l ). Hence, detailed studies of the 
outflow in PKS 1345+12 are important for a number of 
classes of objects (powerful AGN, GPS/CSS radio galax- 
ies and ULIRGs) with the outflow being the best studied 
example of an optical emission line outflow in a ULIRG. 

A key parameter for determining the importance of the 
nuclear outflows in terms of galaxy evolution is the electron 
density of the gas, which allows the calculation of the mass 
outflow rate and kinetic power of the outflow. Traditionally, 
gas densities have been successfully derived for many AGN 
using the [S II]6715,6731 doublet, which requires an accurate 
measurement of the ratio between the two lines in the blend. 
However, sources with highly complex emission line profiles 
can preclude accurate modelling of the [S II] doublet. In the 
case of PKS 1345+12, the velocity widths and shifts of the 
broader components are comparable to the separation of the 
doublet. Hence, whilst total fluxes in each component can be 
established, it is often not possible to determine how much of 
the emission in each component originates from each line in 
the blend. In addition, our previous study of PKS 1345+12 
(H03) has suggested that the densities in the nuclear regions 
of young radio-loud AGN are high (N e > 4xl0 3 cm -3 ), and 
may be high enough to be outside the range of densities to 
which the [S II] diagnostic is sensitive (iV,, = 10 2 -10 3 cm -3 ). 

In this paper, we present new, high quality VLT/FORS2 
spectra of PKS 1345+12, which are both deeper, and cover a 
wider wavelength range, than our previous WHT/ISIS spec- 
tra presented in H03. The new spectra now include further 
density sensitive emission lines, namely the transauroral [S 
II]4068,4076 and [O II]7318, 7319,7330,7331 lines, to pilot a 
new technique to measure the electron density of the circum- 
nuclear emission line gas. Consequently, we calculate for the 
first time the mass outflow rate and kinetic power of the opti- 
cal (warm) outflow in a young radio-loud AGN and ULIRG. 

Throughout this paper we assume the following cosmol- 
ogy: H = 71 km s -1 , D.Q = 0.27 and fi A = 0.73. 



2 OBSERVATIONS AND DATA REDUCTION 

We have obtained new, deep long-slit optical spectroscopic 
observations of PKS 1345+12 in service mode on 22nd 
March 2007. Using FORS2 on the ESO Very Large Telescope 
(VLT) in LSS mode with GRIS-600B+22, GRIS-600RI+19 
and GRIS-600Z+23, and a 1.3 arcsec slit, we obtained spec- 
tra along PA 160 with large spectral coverage (rest frame: 
3200-9200A). 

The data were reduced in the usual way (bias subtrac- 
tion, flat fielding, cosmic ray removal, wavelength calibra- 
tion, flux calibration) using the standard packages in iraf. 
The two-dimensional spectra were also corrected for spa- 
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Figure 1. Rest frame spectrum of the nuclear aperture of PKS 1345+12. The key density diagnostic emission lines used in this paper 
are highlighted. 



tial distortions of the CCD. The final wavelength calibra- 
tion accuracy, calculated uding the standard error on the 
mean deviation of the night sky emission line wavelengths 
from published values (jOsterbrock et all Il996t ) was 0.281, 
0.112 and 0.035 A in the blue, red and z bands respectively. 
The spectral resolution of the night sky emission lines was 
6.65±0.16A, 7.15±0.4lA and 7.02±0.37A in the blue, red 
and z bands respectively. 

Comparison of several spectral photometric standard 
stars taken with a wide slit (5 arcsec) gave a relative flux 
calibration accurate to ±5 per cent. This accuracy was con- 
firmed by comparison to our earlier data (H03). A further 
standard star was observed with a narrow slit, matched to 
that used to observe PKS 1345+12, to correct for atmo- 
spheric absorption features (e.g. A and B bands at ~7600 
and ~6800A respectively). 

The main aperture used was the nuclear aperture - 
1.5 arcsec wide, centred on the nuclear continuum emis- 
sion, and is shown in Figure [T] Before modelling the emis- 
sion lines in the nuclear aperture, the continuum was mod- 
elled and subtracted, a step which has been shown to be 
crucial in the modelling of the broa der, blueshifted compo- 
nents of the emission lines (e.g. H03, iHolt et alJkoOfj , |200St 
2009). Our co ntinuum model com prised a nebular contin- 
uum (see e.g. iDickson et al.l [l995) and a modelled stellar 
contin u um followi n g the meth ods of e.g. iTadhunter et al.1 
l|2005h ; lHoit et al.1 (|2006l . l2007f ). 

The spectra were extracted and analysed using the 
STARLINK packages FIGARO and DIPSO. 



3 RESULTS 

In this paper, our analysis focusses on the density of the 
emission line gas in the circumnuclear regions. We will there- 
fore discuss only the relevant density diagnostic lines ([O II] 
and [S II] ; see Figure [TJ and refer the reader to H03 for a 
detailed discussion of the other strong nuclear emission lines 
in this source. 



3.1 Emission line modelling 

The nuclear emission lines in PKS 1345+12 are strong, 
characterised by highly broadened profiles with strong blue 
wings. A detailed study of the emission line spectrum of the 
nucleus of PKS 1345+12 showed that the emission lines re- 
quired a minimum of 3 Gaussian components to model them 
(H03) . Here, we implement the model derived by H03 for [O 
III]4959,5007 (hereafter the [O III] model). This model com- 
prise^: 

(i) a narrow component, FWHM 340 ± 23 km s -1 at the 
systemic velocity (z = 0.12174 ± 0.00002); 

(ii) a broad component, FWHM = 1255 ± 12 km s" 1 ; 
blueshifted by 402 ± 9 km s _1 with respect to the narrow 
component (note H03 refer to this component as 'interme- 
diate'); 

(iii) a very broad component, FWHM = 1944 ± 65 km 
s _1 , blueshifted by 1980 ± 36 km s _1 with respect to the 
narrow component (note H03 refer to this component as 
'broad'). 

Figure [2] shows the [O III] model applied to the [O 
II]7319,7320,7330,7331 blend, which we have modelled as 
a doublet (i.e. one line centred at ^7320A and one line 
centred at ~7330A), with each line comprising three dis- 
tinct emission line components. The [O III] model gives a 
good fit to the blend. In addition, the [O III] model provides 
reasonable fits to all of the other lines ([S II]6716,6731, [O 
II]3727,3729 and [S II]4068,4076). Note that [O II]3727,3729 
and [S 11)4076,4086 were also modelled as single lines, for a 
number of detailed line-fitting reasons. The measured line 
ratios are plotted on the density diagnostic diagram in Fig- 
ure [3] 

Note that H03 found that the [O III] model did not 
model all of the emission lines perfectly. This included [S 
II]6716,6731 which was best fit by a model with slightly dif- 
ferent FWHM and velocity shifts for the broad and very 
broad components. Furthermore, for all models to [S II], the 

2 The notations narrow, intermediate, broad and very b r oad fo l- 
low the kinematical component definitions of lHolt et al.l {2008). 
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Figure 3. A diagnostic diagram using all of the key density diagnostic lines. The crosses joined by solid lines are the simple, optically 
thick AGN-photoionisation generated using the multi-purpose photoionisation code MAPPINGS. The models are sequences in ionisation 
parameter, U, for optically thick single-slab power-law (F„ <x i/ a ) for spectral index a = -1.0, -1.5, -2.0 and an electron density of n e 
= 100 cm -3 . Also plotted is a density sequence (squares joined by a dot-dashed line) which has a constant power law (a = -1.5) and 
ionisation parameter (U = 0.005) with electron densities of n e = 10 2 , 10 3 , 10 4 , 10 s and 10 s cm -3 . Finally, the squares joined by the 
dashed lines represent the same density sequence but reddened by E(B-V) = 0.5 (faint line) and E(B-V) = 1.0 (bold line) using the 
ISeatonl |l97Sh extinction law. The large red points are the measured (i.e. without reddening correction) data: narrow (triangle), broad 
(square) and very broad (pentagon). In addition, the small red pentagons give the position of the very broad component after applying 
various reddening corrections: A: E(B-V) = 1.44, B: E(B-V) = 0.92 and C: E(B-V) = 0.7. See the text for a discussion of the various 
E(B-V) values. 



intensity ratios for the components needed to be fixed to 
ensure that they were within the permitted range (0.44 < [S 
IIJ6716/6731 < 1.42) and that the fit was physically viable. 
The main reason for the difficulty in modelling [S II] was that 
the widths and shifts of the broader components are compa- 
rable to the separation of the doublet. We find the situation 
unchanged for the new VLT data. For the new diagnostic, 
only the total flux in each kinematic component is required. 
By removing the need for an accurate ratio between the two 
lines in the [S II]6716,6731 doublet, the [O III] model can 
provide a good model to the doublet. To test the robustness 
of our fit, in addition to our 'free intensity' [O III] model 
to the [S 11)6716,6731 blend, we have also tried forcing the 
component ratios so all lie within the range quoted above. 



These give errors of 6%, 9% and 12% on the total fluxes in 
the narrow, broad and very broad components respectively. 

3.2 The electron density 

A key parameter for determining the importance of the nu- 
clear outflows in terms of galaxy evolution is the electron 
density of the gas. As discussed above, due to the highly 
complex emission line profiles in the nuclear regions of PKS 
1345+12, it was not possible to measure the density using 
the [S II]6715,6731 doublet. Here, we pilot a new technique 
based on the transauroral [S II] and [O II] lines. 

Figure [3] shows a new diagnostic diagram con- 
structed for the transauroral density diagnostic 
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Figure 2. The [O II]7319,7320,7330,7331 emission line blend. 
The faint solid line is the observed spectrum, the bold line is the 
overall fit, the dashed lines are the narrow components, the dotted 
lines are the broad components and the dot-dashed lines are the 
very broad components. The vertical line marks the position of 
[Ni III] 7378 in the red wing of the [O II] blend. 

lines [O II]3727/(7318+7319+7330+7331) and [S 
II](4068+4076)/(6716+6731). As discussed above, di- 
agnostics relying only on the total flux in each emission line 
component in each blend bypasses the issue faced with the 
[S II]6716/6731 ratio. In addition, as demonstrated by the 
reddening tracks plotted on Figure [3] whilst the lines span 
a large range of wavelengths, it is not necessary to perform 
separate reddening corrections. The diagnostic is also not 
sensitive to the details of the ionisation model. 

From Figure [3] we obtain electron densities of N e — 
(2.94±?;^) x 10 3 cm" 3 , N e = (1.47lo.l°) X 10 4 cm" 3 and 
N e = (3.16±i'.oi) x 10 s cm -3 for the narrow, broad and very 
broad components respectively. These densities are much 
higher than the original estimates based on [S II]6716/6731 
and suggest that, even without the problems faced mod- 
elling the blend, the diagnostic would not give accurate den- 
sities as the measured densities are beyond the range the [S 
II]6716/6731 ratio is sensitive to. 

Further diagnostic lines sensitive to higher densities 
are [CI III]5518,5534 (sensitive to ~ 10 4 cm" 3 ) and [Ar 
IV]4713,4742 (sensitive to 10 4 - 10 5 cm" 3 ). Both of these 
doublets have been observed in Cygnus A, although they 
are weak (|Tadhunter et al. 1994). Neither doublet was de- 
tected in the original WHT/ISIS spectrum of H03. In these 
data, we detect very faint [Ar IV]4713,4742 although there 
is no evidence of [CI III] emission (see Figure 2J). The detec- 
tion of [Ar IV] further supports the high densities measured 
in PKS 1345+12. In addition, we detect [S III] 9069 (Figure 
1) which also has a high critical density (1.2xl0 6 cm -3 ). 

3.3 Reddening 

As discussed above, PKS 1345+12 is known to harbour 
a dense and dusty circumnuclear cocoon. The level of red- 
dening in the optical spectrum was investigated in detail 
by H03, using three separate techniques: Balmer line ra- 
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Figure 4. The regions of the nuclear rest frame spectra in which 
the higher density sensitive lines [Ar IV]4713,4742A (top panel) 
and[Cl III] 5518, 5534A (bottom panel) would be located. Vertical 
lines mark the wavelengths of the lines. A weak component of [Ar 
IV] 4742 is observed (top panel), but we do not detect any of the 
other lines. 

tios (e.g. Hq/H/3); comparison of the optical results with 
published near-IR data (Paa/H/3) and consideration of the 
nebular continuum as part of the SED modelling. The red- 
dening results of H03 are summarised in Table [1] Here, we 
further investigate the reddening using our density diagnos- 
tic (Figure EJ. 

In addition to electron density sequences, a number 
of red dening (E(B-V)) sequences, derived using the ISeatonl 
l|l979l ) extinction law, are also plotted on Figure [3] From 
Figure |3] the reddening is estimated to be E(B-V) = 0.0 
± 0.1, E(B-V) = 0.7 ± 0.2 and E(B-V) = 0.7 ± 0.2 for 
the narrow, broad and very broad components respectively. 
For the narrow component, this is consistent with the re- 
sults of H03 although the E(B-V) values for the broad and 
very broad components are somewhat higher and lower than 
those derived by H03 respectively. 

In terms of line ratios, a 'cleaner' method would be to 
use unblended line combinations, such as Paa/H/?. H03 at- 
tempted to use Paa/H/3 to improve the reddening estimate 
but issues such as differing components and uncertainties 
about possible slit losses between the optical and near-IR 
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Table 1. Reddening values and emission line data. Columns are: (a) Quantity, where: 
Ha/H/3 is the flux ratio between the measured values of Ho and H/3; E(B-V)(Ha/H/3) 
is the E(B-V) value calcula t ed us ing the Hcr/H/3 ratio and the standard interstellar ex- 
tinction curve from [SeatorJ l|l979f) ; Fjjff is the reddening corrected H/3 line flux (10 -15 
erg s _1 cm -2 ); L H ^ is the rest frame H/3 luminosity (10 40 erg s _1 ); E(B-V) (limits) gives 
the lower limits on the E(B-V) value derived for the broad and very broad components 
by H03; E(B-V) (diagnostic) is the E(B-V) value derived from the diagnostic diagram in 
Figure [3] (ft) value of particular quantity for the narrow component; (c) uncertainty for 
narrow component; column pairs (d) and (e) and (/) and (g) are as for (ft) and (c) but 
for the broad and very broad components respectively. 

a derived from considerations of the nebular continuum (H03). 6 derived from Ho/H/3 
using the model to the Ha,[N II] blend which gives the lowest possible flux in Ha (H03). 
c derived using the E(B-V) values given in row 2 of this Table. d luminosity from H03 
updated to the cosmology used in this paper. e derived using the new (VLT) H/3 fluxes 
and E(B-V) (diagnostic) values. 
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data proved difficult. We have since obtained the unpub- 
lished near-IR spectrum of Veilleux (priv. comm.) in an 
attempt to model Paa with our three component [O III] 
model. Whilst we can obtain a good fit, it is no better than 
the original two component model, and still does not help 
with regards to the issue of slit losses. 

Given the uncertainties regarding the degree of redden- 
ing, H03 attempted to derive limits for the reddening using 
a more indirect technique. By generating nebular continua 
with various E(B-V) values, H03 set a lower limit on the red- 
dening in the broad component, the component which dom- 
inates the flux in the nebular continuum. This technique 
gave a lower limit of E(B-V) > 0.3 for the broad compo- 
nent, which is consistent with all reddening estimates for 
this component. 

Given the uncertainties in the reddening estimates, we 
have also plotted the reddening corrected line ratios for the 
very broad component using various E(B-V) values, in or- 
der to assess how the degree of reddening affects our derived 
density. For this exercise we have chosen three E(B-V) val- 
ues: 

(i) point A: E(B-V) = 1.44; the value derived by H03 
using Ha/H/3; 

(ii) point B: E(B-V) = 0.92; the lower limit for the red- 
dening derived by H03; 

(iii) point C: E(B-V) = 0.7; the value derived by our di- 
agnostic diagram. 

As can be seen from Figure [3] whilst the poorly constrained 



E(B-V) value will affect the reddening corrected H/3 lumi- 
nosity, it has little effect on the derived electron density. We 
therefore conclude that the electron densities derived using 
this technique are robust. 



4 DISCUSSION 

4.1 How important is the warm outflow in PKS 
1345+12? 

In recent years AGN feedback, in the form of quasar-driven 
outflows, has been relied upon to explain a number of 
observed relations, such as the close correlation between 



Silk & Reesj 19981; Fabian 


19991; Ferrarese & MerrittJ 


200d; iGebhardt et all 12000 


; iTremaine et al.l 


2003; 
2005; 


Marconi & HuntJ 120031; iDi Matteo et all 


Kurosawa et all I200SL Booth & Schavd 2009) and the 



high end shape o f the galaxy luminosity function (e.g. 
iBenson et al. 2003). Currently, the majority of AGN feed- 
back models require a significant amount of the available 
accretion energy to power the outflow s (~5-10%; e.g. 
| Di Matteo et alj|2005l; ISpringel et al.ll2005l ; iBooth fc Schawl 
2009; Kurosawa et alj|20090 although this has recently been 
brought into question, with much lower energy injection 
required if a t wo-stage feedback mo del is implemented 
(~0.5% of Lhoi; lHopkins fc Elvidl20ich . 

PKS 1345+12 provides an ideal test case for investi- 
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gating the impact of AGN feedback in the early stages of 
AGN evolution. Clearly currently still undergoing a major 
merger/interaction (c.f. double nucleus within a single dif- 
fuse halo with a number of tidal features and evidence for 
recent star formation), the small, recently triggered radio 
jets are currently expanding through a dense natal cocoon 
of gas and dust (c.f. large electron densities and reddening, 
high X-ray and HI absorbing columns). Our optical spec- 
tra reveal clear signatures of AGN driven outflows in the 
strong, highly broadened and blueshifted optical emission 
lines. Whilst the outflows are fast, it is important to com- 
bine this kinematic information with other properties, such 
as the gas density, to determine the key parameters for un- 
derstanding the importance of AGN feedback - the mass 

outflow rate and the kinetic po wer. 

Following the derivation in iHolt et al l (|2006l ). the mass 
outflow rate for a spherical outflow can be expressed as: 

M — 3L(Hf3)rn p v ou t (1) 

Nof^hVHfjT 

where L(H/3) is the H/3 luminosity, m v is the mass of the 
proton, Vout the outflow velocity, N is the electron density 
of the ga s, a e Jg is t he effective H/3 recombination coeffi- 
cient (see lOsterbrocklfl 9891 ). Vup the frequency of H/3, h the 
Planck constant and r is the radius of the spherical volume. 

We take H/3 luminosities measured from our new VLT 
spectra, modelled with the [O III] model of H03, and cor- 
rected for extinction using the E(B-V) val ues deri v ed by 
our diagnostic diagram (Figure [3J and the ISeatonl (|l979l ) 
extinction law, all of which are are summarised in Table [T] 
In our currently adopted cosmology, these are: L(H/3)n = 
(4.3±0.7)xl0 40 erg s" 1 ; L(H/3) B = (1.55±0.18) x 10 42 erg 
s _1 and L(H/3)vb = (5.6±0.15) x 10 41 erg s" 1 . The rest 
frame outflow velocities (v out ) axe vb = 402±9 km s _1 
and «vb = 1980±65 km s~ . Using our new diagnostic 
tool (Figure [3|, we measure electron densities of N e = 
(2.941?;^) x 10 3 cm" 3 , N e = (1.47^'j?) x 10 4 cm" 3 and 
N e = (3.16ti.o?) x 10 s cm -3 for the narrow, broad and very 
broad components respectively. If we assume case B recom- 
bination theory for an electron temperature T = 10, 000 K 
and that the outflow is spherical, with a radius half of the 
diame ter of the [O HI] emissi on line region imaged by HST 
ACS (Bat cheldor et al.1 |2007r i (i.e. 7.5 milli-arcseconds or 
r = 162pc for our cosmology), we calculate mass outflow 
rates of M = 7.1±f% M yr _1 and M = 0.6±gJ M yr _1 
in the two outflowing components respectively, with a to- 
tal mass outflow rate of M = 7.Ttl\l M y r_1 - Tnis es " 
timate is comparable to the range of mass outflow rates 
calculated for another young radio source, PKS 1549-79 
(H06) but lies between the mass outflow rate s estimated for 
the N LR of Seyferts (0.1-1.0 M© yr" 1 ; e.g. IVeilleux et ail 
120051 ) and the much larger neutral outflows detected in 
other ULIRGs (10 2 -10 3 M© yr _1 ; e.g. lRupke et al.ll2005al lbh. 
We estimate that the total mass in the warm gas out- 
flow is Moutflow = (Stfj) x 10 5 M with filling factors^ of 
e = (4.4ti;s) x 10 ~ 4 and 6 = 0- 6 -o.l) x 10 ~ 7 for the regions 
emitting the broad and very broad components respectively. 

In order to gauge the impact of the warm gas outflow 
on the circumnuclear gas, it is important to estimate the 

3 See equation 3 in lHolt et all |2008). 



kinetic power of the outflow, including both the radial and 
turbulent components in the gas motion. Assuming that the 
relatively large linewidth of the outflowing gas reflects a tur- 
bulent motion that is present at all locations in the outflow 
region, the kinetic power is: 

M 

E = 6.34 x 10 35 y(?w + (FWHM) 2 /1.85) ergs -1 (2) 

where FWHM is the full width at half maximum of the [OIII] 
line component in km s , and M is the mass outflow rate 
expressed in solar masses per year. The kinetic powers are 
E = (2.3±0.9)xl0 42 ergs" 1 and E = x 10 42 ergs" 1 

for the broad and very broad components. Hence, the total 
kinetic power in the warm gas outflow is E = (3.4i];"|) x 
10 42 erg s" 1 . 

Combining the kinetic power with the bolometric lu- 
minosity provides an estimate of the fraction of the total 
available power fr om accretion onto t he black hole driv- 
ing the outflows. iDasvra et al.l |2006) publish two black 
hole mass estimates for the active western nucleus of PKS 
1345 +12, one derived us ing M B h - o as M B h = 6.54xl0 7 
M© (|Dasvra et alj I2OO6I . and a smaller mass estimate of 
Mbh = 1.75xl0 7 Mq based on the assumption that |Lfir 
= Lboi = LEdcEl- Comparing the two black hole masses gives 
an idea of the Eddington ratio: 

Mbh, fir _ -fcboi _ n 97 ,„\ 
~M — 7 — ^ > 

A'BH.dyn -t^Edd 

This estimate for the Eddingt on ratio is reason able. For the 
similar source PKS 1549-79, IHolt et all (120061 ) derived Ed- 
dington ratios in the range 0.3 < Lboi < 5 7 depe nding on 
the black hole mass used. Furthermore, IWul l|2009l ) collated 
data from the literature for a large sample of compact radio 
sources. Of the 51 sources with quoted Lboi/LEdJj, the Ed- 
dington ratio is typically high (10 -2 < Lboi/LEdd < several) 
compared to samples of radio loud and radio quiet quasars 
and broad line radio galaxies (typically Lboi/ LEdd < 0.1; 
e.g. iDunlop et alj|200a iMarchesini et nl.ll200 lls but similar 
to th e ratios observed in narrow line seyfert Is (Lboi/LEdd = 
0.16: lBian et alfeOOcf ); with L b oi/L E dd = 0.28 and L bo i/L E dd 
= 0.26 for the CSS and GPS sources respectively. The es- 
timate of Lboi/L E dd = 0.27 for PKS 1345+12 derived from 
the ratio of the two black hole estimates is therefore entirely 
consistent with other GPS sources where this ratio is known 
with more certainty. 

With this in mind, we calculate the fraction of the accre- 
tion power in the outflow as: E/L\,o\ = (1-3 ± 0.2) x 10~ 3 . 
Hence, if the assumed bolometric luminosity conversion is 
correct, only a small fraction of the available accretion 
power is driving the warm gas outflow. This is similar to 



4 It is also possible to estimate the bolometric luminosity from 
the extinction-corrected [O HI] emission line lum inosity using the 
conversion factors given in Lamastra ct al. Our reddening 
corrected total [O III] luminosity is L [onl ] = (2.11 ± 0.25) x 10 43 
erg s — 1 and corresponds to a bolometric luminosity of L[oiii] = 
(9.5 ± 1.1) X 10 45 erg s" 1 , which is similar to the estimate derived 
from the FIR luminosity. 

5 Whilst PKS 1345+12 is included in the paper, IWul 1120091) di d 
not include the estimate of Lboi/LEdd from lDasvra et al.l (2006). 
which is expanded for clarity above. 
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the findings for anot her compact radio source, PKS 1549- 
79 (|Holt et alj 2006). This is in stark contrast to the re- 
quirements of the majority of quasar feedback models which 
require a much larger fraction of the accretion power of 
the black hole (~5-10% of Lbol) to be thermally coupled 
to the circumnuc l ear g as, but similar to the recent work of 
iHopkins fc Ehdsl (|2010h . 

It is also unlikely that the warm gas outflow is ca- 
pable of removing all the warm/cool gas from the cen- 
tral regions of the host galaxy by itself. If we assume that 
PKS1345+12 has a total mass Af to tai ~ 10 11 M Q , and 
a gas mass M gas ~ 10 10 Mq contained within a radius 
of 5 kpc - conservative assumptions for ULIRGs which 
are likely to have a larger gas mass concentrated within 
a smaller radius - the gravitational binding energy of the 
gas is Ehmd ~ GM gas Mtotai/-Rgas ~ 2 x 10 58 erg. In com- 
parison, the warm gas outflow will deposit only ~ 10 56 erg 
into the surrounding ISM, assuming that it can persist in its 
current form for the typical 10 7 yr lifetime of an extragalac- 
tic radio source. Again, these observations are consistent 
with the limits calculated for another compact (young) ra- 
dio source PKS 1549-79 (HOB). This result may suggest that 
PKS1345+12 will grow into a large radio galaxy in a gas-rich 
host. However, with its double nucleus, PKS 1345+12 is still 
in the pre-coalescence stage of the merger, observed during 
the first peri-centre passage, but immediat ely before (within 
~0.1 Gyr) of the coalescence of the nuclei (|Tadhunter et all 
2010). According to the merger simulations, PKS 1345+12 is 
therefore yet to pass through the most active phase (both in 
terms of nuclear activity and star formation) of the merger 
when the nuclei coalesce. An increase in nuclear activity may 
increase the strength of the winds/outflows and an increase 
in star formation will also consume large quantities of the 
ambient ISM. It is therefore unclear whether or not PKS 
1345+12 is a progenitor of a gas rich radio galaxy. 

It is clear that the kinetic power of the warm gas out- 
flow in PKS 1345+12 is significantly less than that required 
by the majority of current quasar feedback models. Could 
this be because we have not accurately measured all of the 
input parameters or made inaccurate assumptions? In the 
above analysis, the key observable parameters are the elec- 
tron density, the size of the outflow region, the luminosity of 
the H/3 line and the black hole parameters (black hole mass 
and the Eddington ratio). Following through, an underes- 
timate of the kinetic power, E, would result from an over- 
estimated electron density and/or outflowing region size, or 
from an underestimated H/3 luminosity. Furthermore, if the 
black hole mass and/or the bolometric luminosity is overes- 
timated, the impact of the outflow may be underestimated. 

We are confident that the size of the outflowing region 
has not been overestimated as recent high-resolution HST 
imaging of PKS 1345+12 have shown the optical line emis- 
sion to be on a similar scale to the compact radio emission 
ijBatcheldor et alJl2007l ). If anything, the outflowing region 
may be significantly larger if projection effects are impor- 
tant. 

The electron density derived by this new technique is 
high (N e ~ 10 4 " 5 cm" 3 ) and therefore one to two orders 
of magnitude above typical NLR densities in the literature 
measured using traditional diagnostics. As discussed in Sec- 
tion 3.3, the large uncertainties in the reddening estimate 
will have a significant effect on the line ratios. However, as 



demonstrated by Figure [3l even with large changes in the 
reddening for the very broad component, the derived den- 
sity remains unchanged. Despite this, it is still important to 
consider whether such high densities are realistic and what 
effect this has on our overall conclusions. 

Whilst it is possible the new technique overestimates 
the density, very high densities are not surprising; high gas 
densities ha ve been observe d in a number of compact ra- 
dio sources (|Holt et al.ll2009T ) and the original study of PKS 
1345+12 using the [S II]6716/6731 ratio calculated a rela- 
tively high lower limit (jV e > 5 x 10 3 cm -3 ; H03). Further- 
more, large quantities of gas in the nuclear regions of PKS 
1345+12 have been observed in a variety of wavebands (see 
Section 1). Using the lower limit for the electron density 
provides upper limits for the mass outflow rate (M < 21 
M y r_1 )i tne kinetic power (E < 3.8 x 10 43 erg s" 1 ) and 
the fraction of the available accretion power in the outflow 
{E/L Edd < 4 x 10" 3 ). Hence, even if our electron density is 
overestimated by as much as two orders of magnitude, the 
amount of energy in the outflow is still significantly lower 
than the required 5-10%. 

It is possible that the errors on the derived H/3 luminos- 
ity are large. As discussed above, the nuclear regions of PKS 
1345+12 are known to contain large quantities of gas and 
dust. The optical emission line components are known to be 
heavily extinguished (H03 and Section 3.3), whilst radio and 
X-ray studies detect large absorbin g columns: Nu > 10 22 
cm" 2 (radio; iMorganti etail 120051 ) and A H = 2.3 x 10 22 
cm" 2 (X-ray; Siemiginowska priv. comm.). However, as 
demonstrated in Section 3.3, it is difficult to constrain the 
degree of reddening, although all estimates discussed in Sec- 
tion 3.3 do agree within the (large) errros. Above, we cal- 
culated the various outflow properties, such as mass outflow 
rates and kinetic powers, based on the H/3 luminosities de- 
rived from the new H/3 flux measurements from the VLT 
spectra, and the reddening estimates derived from the di- 
agnostic diagram in Figure [3] We now calculate the same 
quantities using the new H/3 fluxes and the original redden- 
ing estimates from H03 in order to understand the possible 
uncertainties in our results. 

In short, for this case, we calculate the following: (i) 
mass outflow rates of M B = 2.7 M© yr" 1 and M VB = 7.2 
M© yr" 1 with a total mass outflow rate of Mtotaj = 9.9 
MQ yr" 1 ; (ii) volume filling factors of ee = 1.7xl0~ 4 and 
£vb = 1.9xl0" 6 ; (iii) total gas mass of Mtotal = 8xl0 5 
MO; (iv) kinetic powers of E B = 8.7 x 10 41 erg s" 1 and 
_Evb = 1.4 x 10 43 erg s" 1 with a total kinetic power of 
Etotsi = 1.5 x 10 43 erg s" 1 ; (v) fraction of the total avail- 
able accretion in the warm outflow: -E/LEdd = 6xl0~ 3 . 
Hence, implementing the original reddening estimate from 
H03 rather than those derived from Figure [3] decreases the 
impact of the broad component by a factor of ~3, but boosts 
the impact of the very broad component by a factor of ~12. 
Overall, this increases the total impact of the outflow by a 
factor of 4.6, from 0.13% of Lb ; to 0.6% of Lb ; and demon- 
strates that our results are robust against large variations 
in the reddening. 

The final possible source of uncertainties relates to the 
black hole properties, namely the mass, the bolometric and 
Eddington luminosities and the Eddi ngton ratio . In the 
above calculations, we have followed iDasvra et ail (|2006l ) 
and assumed that tjLfir = Lboi = LEdd provides a good 



Impact of the warm outflow in PKS 1345+12 9 



estimate of the bolometric luminosity. Whilst this conver- 
sion may be a good general assumption, in highly complex 
sources like PKS 1345+12, the FIR emission is likely to con- 
tain many com ponents. We have a lso compared the black 
hole masses in iDasvra et alj (|2006T ) and estimated an Ed- 
dington ratio of 0.27. To quantify the possible errors on 
our result, we take the extreme case and assume that PKS 
1345+12 is accreting at the Eddington ratio, with the larger 
of the two black hole masses. In this case, we calculate 
E/L Edd = (3.6+0.16) x 10~ 4 . Hence, potential errors in 
the black hole mass are likely to reduce the impact of the 
outflow. 

Given t he above, a s sugge sted for a similar source PKS 
1549-79 by iHolt et alj {2006), it is likely that the warm 
gas outflow detected in the optical emission lines accounts 
for only a small fraction of the total outflow with the re- 
maini ng mass in ho t ter or cooler phases of the ISM. In- 
deed iMorganti et alj (|2005T ) reported the discovery of mas- 
sive neutral outflows in a sample of 7 compact radio sources 
which includes PKS 1345+12. The HI outflow in PKS 
1345+12 is: M ~ 8 - 21 M yr"\ with the lower end of 
the range similar to that for the warm gas outflow. Taking 
this range of mass outflow rates, and assuming the outflow 
velocity and FWHM to be 600 km s" 1 (i.e. FWZI/2; Mor- 
ganti priv. comm.), we find the HI outflow has a kinetic 
power of 1 x 10 42 < E (erg s" 1 ) < 4 x 10 42 which accounts 
for 4 x 10~ 4 < E/L hol <lx 10 -3 of the available accretion 
energy. 

Combining the HI and optical results, the total observed 
outflow is currently of the order of 16-29 Mq yr _ , with 
a kinetic power of 4.4xl0 42 < E < 7.4xl0 42 erg s _1 and 
accounts for 0.2-0.3% of the available accretion energy, Lbol- 
It is also interesting to consider the amount of energy 
locked into the radio emission by comparing the radio power, 
Qjet, to th e accretion energy. Taking the radio jet power cal- 
culated by|Wu| Hil) of log(Q ict ) = 44.6fl and the above 
bolometric luminosity, we find Qjet/^boj = 0.16. In other 
words, around 16 per cent of the total available accretion 
power is currently powering the radio jets. Whilst there is 
various evidence to suggest that the HI and optical emission 
line outflows are driven by the r adio jets (e.g. high out- 
flow velocities e.g. IHolt et al"1l2008l ; co-spatial radio and op- 
tical emission on similar scales in GPS and CSS sources and 
in some CSS sources, evidence for radio-optical alignment 
Ide Vries et alJll999l:|Axon et all200d;IOT)ea et al.ll2002l; 



Labiano et al. I 120051 ; Batchcldor et al. I 120071 ; iPrivon et al.1 
20081 ). it appears that the mechanical coupling of the radio 



jet to the ISM is relatively inefficient. 

PKS 1345+12 is also detected in X-rays (e.g. 
lO'Dea et af]|2000l ; Isiemiginowska et alj|2008h with a soft X- 
ray excess (L(0.5-2keV) = 2xl0 42 ergs sec -1 ) and a harder 
X-ray spectrum than typical for AGN (Siemiginowska priv. 
com.). The emission has two clear features: i) an extended 
(~15 kpc) diffuse component, which may be the result of 
thermal emission in the host galaxy, and ii) an elongated 
feature to the SW of the active nucleus, which appears to 
be closely related to the VLBI radio source. Further obser- 
vations are required to establish whether a hot gas outflow 
is present in PKS 1345+12. 



6 At 178MHz. 



To summarise, after consideration of all possible sources 
of error, the most extreme kinetic power we calculate is still 
at least an order of magnitude smaller than that currently 
required by the majority of galaxy evolution models but 
similar to the energy requiremen ts suggested by recent w ork 
on a two-phase feedback model (|Hopkins fc Eividlanoh . 



4.2 PKS 1345+12 — a young radio source? 

It is now well accepted that compact (GPS and CSS) ra- 
dio sources are compact due to evolutionary stage rather 
than due to the effects of frustration by a dense ISM (e.g. 
iFanti et al.1ll995h . However, recent simulations of radio jet- 
ISM interactions have shown that, if the ISM forms dense 
clouds, much smaller masses could significantly distrupt 
the propagation of the jet. For example, jet interactions 



with dense clouds (e.g. of the ord er M a 



10° M B ) 



may be significantly disrupted (e.g. ISutherland fc Bicknelll 
2007) compared to the much larger gas masses required 
(M gas ~ 10 9 " 11 M ) required if the gas is distributed in 
a homogeneous or clu mpy distribution fe.g. lDe Young|[l993l ; 
ICarvalho|[l99l . 1 19981 ) . 

We have seen that the nuclear regions of PKS 1345+12 
harbour a rich, dense ISM, but is there sufficient mass to 
confine and frustrate the radio source? 

In our original study in 2003, using the following equa- 
tion: 



M„ 



L(H/3) 



(4) 



we derived limits on the gas mass for PKS 1345+12 of: 
M gas > 2.61 xlO 5 M Q , Mgas < 0.92xl0 5 M and M gas < 
5.64x10 Mq in the regions emitting the narrow, broad and 
very broad components respectively; with an upper limit of 
order 10 6 M for the total mass of line emitting gas in the 
kinematically disturbed emission line components. 

Using the new gas densities and H/3 luminosities mea- 
sured in this paper, we now derive gas masses of: M gas = 
(1.2t^) xlO 5 M , M gas = (7.7l| 1 )xl0 5 M and M gas = 
(2.9±i'°)xl0 4 M Q for the narrow, broad and very broad 
components respectively. This gives a total gas mass of 
(9.2±3.o) xlO 5 M Q with (8+3) xlO 5 M in the outflowing 
components. 

H03 discussed that, even with the original estimate of 
<10 6 M Q , it was unlikely that there was sufficient warm 
gas in the nuclear regions of PKS 1345+12 to confine and 
frustrate the young, expanding radio jets and only the pres- 
ence of dense clouds would be able to significantly disrupt 
the propagation of the jet. The newly derived gas densi- 
ties suggest the warm gas mass is one-two orders of mag- 
nitude smaller than the previously calculated limit. It is 
therefore unlikely that the propagation of the radio jets 
in PKS 1345+12 will be significantly affected by the warm 
ISM in this source. However, it should be noted that much 
higher gas mas ses have been meas ured in e.g. molecular gas 
(3.3xl0 10 Mfl: lEvans et al.lll999l ). 



5 CONCLUSIONS AND FUTURE WORK 

Using new, deep VLT/FORS spectra of the young radio 
source PKS 1345+12, we have piloted a new technique using 
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the transauroral [S II] and [O II] emission lines to determine 
the electron density in a source where the highly broadened 
and complex emission line ratios precluded measurements 
using the traditional diagnostics (e.g. [S II]6716/6731). Our 
observations reveal that the outflowing warm gas in this 
compact radio source is dense {N e ~ 10 4 " 5 cm' 3 ). In 
addition to previous issues regarding the modelling of [S 
II]6716,6731 doublet, this shows that the electron densities 
in this source are much higher than the range of densities to 
which the diagnostic is sensitive. 

The newly determined electron densities have enabled 
us to calculate key parameters for determining the impact 
of the fast, warm outflow in this source (mass outflow rate, 
kinetic power). Whilst at face value the warm gas outflow 
appears extreme, with outflow velocities of up to 2000 km 
s _1 , the outflow is driven by only a small fraction of the 
energy available from accretion power. Comparisons with 
the majority of AGN feedback models in the literature sug- 
gest that the outflow in PKS 1345+12 is at least an order 
of magnitude smalle r than required, However, recent work 
by I Hopkins fc Elvis! (|2010l ) suggests that with a two-stage 
feedback model, the initial feedback requirements may be 
similar to what is observed in PKS 1345+12. 
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